Decoherence of mesoscopic states of cavity fields 
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Abstract 



We show that two-atom correlation measurements of the type involved in a recent ex- 
perimental study of the evolution of a mesoscopic superposition state prepared in a definite 
mode of a high-Q cavity can be used to determine the eigenvalues of the reduced density 
matrix of the field, provided the assumed dynamical conditions are actually fulfilled to ex- 
perimental accuracy. These conditions involve i) a purely dispersive coupling of the field to 
the Rydberg atoms used to manipulate and to monitor the cavity field, and ii) the effective 
absence of correlations in the ground state of the system consisting of the cavity coupled to 
the "reservoir" which accounts for the decoherence and damping processes. A microscopic 
calculation at zero temperature is performed and compared to master equation results. 
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The conception of an experimental setup to produce a superposition of mesoscopic quan- 
tum states and to observe its subsequent decoherence |IJ in microwave cavities has been pre- 
ceded by a series of meticulous experimental and theoretical studies 10 -[fj. Several schemes 
have been proposed, based on two atom correlation measurements involving circular Rydberg 
atoms, and the actual implementation of one such scheme has been recently reported It 
can be broadly pictured as involving two successive stages, which consist respectively in the 
preparation of the desired superposition state and in the subsequent monitoring of its time 
evolution under influence of the unavoidable loss processes present in the cavity. As the last 
step in the preparation stage, what is now often referred to as the "reduction postulate" 
is used instrumentally in an essential way to select the desired state of the field (see below). 
The subsequent stage, however, involves just strictly causal quantum evolution laws, and is 
therefore, in principle, amenable to a complete theoretical description, including the coher- 
ence loss of the prepared superposition state @, ||. In particular, as it is appropriate for 
bona-fide physical properties and processes, the coherence loss should be characterized and 
described in a basis independent fashion. We show in continuation that this can be done in 
a particularly simple way in the present case, given the specific properties of the prepared 
superposition state, and no more than ordinary assumptions concerning the external cou- 
plings underlying the mechanism of coherence loss. More specifically, we show that, under 
the experimental conditions stated in || , two-atom correlations of the type involved in the 
measurements reported there can be used in order to determine completely the eigenvalues 
of reduced density matrices which describe the possible quantum states of the cavity field 
following the preparation stage. The time dependence of the observed two-atom correlations 
subsequently to the preparation of the superposition state constitutes therefore a measure- 
ment of intrinsic features of the evolution of the field which bear directly on the coherence 
loss process. Moreover, to the extent that this evolution depends both on the nature of the 
initial state and on dynamical assumptions concerning the dissipative process, through an 
analysis of the type proposed here such measurements can provide experimental checks on 
the actual outcome of the preparation stage and on the adequacy of the assumed dissipative 
dynamics, at no higher cost than a less inclusive analysis of the two-atom correlation data. 

The experiment reported in || involves a high-Q cavity C located between two low-Q 
cavities (Ramsey zones Ri and R2) fed with classical fields which steer the internal state of 
the crossing circular Rydberg atoms. The transition between two near atomic levels, usually 
denoted as | e) and | g), is resonant with the fields in cavities R\ and R2. Each atom is 
initially prepared in the state | e). After leaving R\ it is in a superposition of the states 
I e) and | g). The fields in Ri and in R 2 (and their relative phases) are chosen so that their 
action on the atoms is given in both cases by the unitary transformation 

Ur: |e>--^(|e>+|0»; I 9) - ^= (- I e>+ | g)) . 

The high-Q cavity C stores a coherent field | a). The coupling between the atom and the field 
in this cavity is measured by the Rabi frequency Q. Two situations will be considered: a) 
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the frequency u of a relevant mode of C is nearly resonant with the transition frequency U0i e 
connecting | e) to a third state | i) (assumed for defmiteness to lie above | e) in energy) and 
is far off resonance with any transitions involving level | p) 0; and b) u is nearly resonant 
with the transition frequency u eg connecting | e) to | g) 0, |6|. In both cases the detuning S 
is assumed large enough so that real transitions are effectively hindered and the atom-field 
interaction leads essentially to 1/5 dispersive frequency shifts. In this way, the atom- field 
coupling during a time t produces an atomic level dependent dephasing of the field and thus 
generates entangled states of the atom and field. This can be represented by the unitary 
operators 

case a : = e~ l ^ a \ e) (e | + | g) (g \ 

(1) 

case b: = e i(l){ata+1) \ e) (e \ +e~ l ^ a \ g) (g \ 

with cf)(t) = Q 2 t/5. In these expressions the operators a\ a create and annihilate photons 
in the relevant mode of C. The effect of these operators on a coherent state of this cavity 
amounts thus simply to a change of phase which depends on the state of the intervening 
two-level atom, so that the resulting state of the field can be expressed in terms of just two 
coherent states. It should be emphasized that this feature depends crucially on the validity 
of the purely dispersive limit mentioned above, which requires values of Q/5 sufficiently small 
to ensure the validity of the implied low-order perturbative treatment [H]. 

This idealization of the involved dynamics allows one to describe the experiment in a 
symbolically compact form. When an atom is sent through the apparatus, the change in the 
state of the combined system consisting of atom plus field is described by UrUcUr, so that 
for a factorized initial state p F ® \ e)(e | (with p F —\ a) (a |) one has 

U R U c U R (p F ® |e)(e|)444 

for the combined system. The field state resulting from the final observation of the atom in 
state | e) or | g) is obtained by replacing the product of unitary operators respectively by 
the reduced operators 



U % = ( e g \ U R U C U R I e) (2) 

(or their adjoints) with an appropriate trace normalization factor. This is what we referred 
to as the preparation stage of the experiment. The instrumental role attributed to the 
"reduction postulate" in this stage appears explicitly in this last step. One gets, for the two 
cases considered, 



U 



(a) 



=Fl 



and U, 



(b) 



1 r 
2 



(3) 
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The indices e, g and the double signs refer respectively to the final observation of the atom 
in states | e) and | g). 

We now turn to the second stage of the two atom correlation measurement. The second 
atom is sent through the system following the detection of the first. The passage of this 
second atom trough the apparatus and its detection can be described using the same formal 
tools adopted for the first atom. Let pf be the trace normalized reduced densities which 
describe the state of the field in C after the detection of the first atom in e or g. The 
conditional probabilities for the various possible results obtained in the final measurement 
of the state of the second atom, for given outcome of the measurement of the final state of 
the first atom, are respectively given by 



Pee = Tr 



UlU e p F ]; P nP = Tr\ulU, 



F 



P eg = Tr [up g p F ] ; P gg = Tr \u%p 



(4) 



Note that in e.g. P ge the first atom is detected in state g and the second in state e. Because 
of the assumed normalization one has P ee + P eg = 1 = P ge + P gg , and the correlation signal 
measured in ref. H can be written simply as 



r, = P ee - P ge = Tr \u\U e (pf - p F g )\ . (5) 

The reduced densities appearing here refer to the state of the field in C at the time of the 
passage of the second atom through the apparatus. On the one hand, these expressions 
connect the two-atom correlation data directly with theoretical models for the fully causal 
time evolution of p F during the time lapse between the passage of the two atoms. This 
involves in fact the only truly "complex" aspect of the dynamics, which is responsible for 
the dissipative evolution of the field state in C. Since moreover one can always express the 
reduced densities in terms of their (time dependent) eigenvectors and eigenvalues, Eqs. (|) 
and (|J) also show directly that measured values of the various conditional probabilities and 
of 1] can be expressed in terms of these same ingredients. As shown below, the eigenvalues 
are in fact determined by P ee and P ge , provided specific but ordinary assumptions concerning 
the dissipation mechanism are made. In this sense, they provide a basis independent char- 
acterization of the decoherence process much as quantities based on more general constructs 
which are nonlinear in p such as e.g. the "idempotency defect" 1 — Trp 2 or entropy-like 
functionals of the density f7j. The assumptions concerning the damping mechanism will 
allow, in particular, for an exact treatment, independently of the also customary recourse to 
master equation approximations. 

In order to model the dissipative effects in C during the time between the passage of 
the two atoms, we introduce as usual a set of harmonic oscillators with frequencies Uk and 
creation/annihilation operators written as b k , (the "heat bath") coupled to the relevant 
mode of the cavity through the bilinear expression 
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Hint = ^1k{a ] h + b\a). (6) 

k 

This differs from the also often used position-position coupling by the omission of "anti- 
resonant" terms involving products of two creation and of two annihilation operators. These 
terms are actually expected to contribute very little to the damping process in view of the 
smallness of the ratio 7^/0; in the present context. The resulting system of coupled oscillators 
has the well known properties that i) its ground state is the product of the ground states 
of the different oscillators (no ground state correlations) and ii) an initial state written as a 
product of coherent states, i.e. 

|t = 0>=|a>ni Pk) (7) 
k 

retains this form at all times, the amplitudes a(t) and /?&(£) being solutions to the linear 
coupled equations 

ia = J2 7*A 

k 

(8) 

ifi k = (u k - u)(3 k + 7 fc a 

with the given initial conditions, in an interaction picture where the time dependence asso- 
ciated with the constant of motion Uuj {^a + J2k b\bkj has been eliminated. This solution 
is in fact the relevant one when we then take as initial states the states generated by the 
action of the operators obtained in Eq. (§) on a coherent state | a) and a heat bath at 
temperature zero. Since the state of the relevant mode of C is then the sum of two lin- 
early independent (though, strictly speaking, never orthogonal) coherent states, we see that 
we will have at all times a sum of two terms, each being of the type shown in Eq. ([?]). 
Consequently, the reduced density describing the state of the field in C will have its trace 
exhausted in a (time dependent) two-dimensional subspace of the corresponding oscillator 
space This particular feature is in fact retained in the currently used master equation 
description of the damping process (see e.g. ||), which however gives quantitatively differ- 
ent results both at very short times (where it gives a linear rather than a quadratic time 
dependence for decoherence measures such as the idempotency defect 0) and for long times 
in the asymptotically damped region. 

The computationally simplest case to consider corresponds to case a) as defined above, 
with <p = 7r. The "initial" state (after the detection of the first atom) is then 

\t = 0) = N a [\ a)=F I -a)} ® | 6 ) 
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where | Ob) stands for the "bath" ground state (zero temperature), the double sign corre- 
sponds to the detection of the first atom in e or g and N a is the appropriate normalization 
factor. The solution of Eqs. (0) then gives at time t 



t) = N Q 



a(*)>III &(0>T|-a(t)>ni-&(*)> 

k k 



It is then a simple matter to obtain the C-field reduced density at time t as 



Pf(0 = K [| a(f)><a(f) | + | -a(t)){-a{t) | T r 6 (t) (| a(f))<-a(f) | + | -a(t))(a(t) |)] (9) 



where the real function r&(t), which contains the decoherence effects of the bath oscillators, 
is given by 

r b (t) = H(-i3 k (t)\f3 k (t)). 

k 

The non- vanishing eigenvalues \±(t) of each of these two matrices are 

,<$>,,, (i + r„(t)) (i t r t (t)) (j) (i-r a (t))(i±r t (t)) 

A+(t)= 2(i±r.(0)) ; A - (i,= 2(i±r„ ( o)) (10) 

where T a (t) = (—a(t) \ a(t)), with associated eigenvectors 

| A±(t)) oc | a(t))± | -a(t)). 

Since in this special case the operator exp(i7ra^a) appearing in U.\U e is simultaneously diag- 
onal (with eigenvalues ±1 for | A±(t))) with the reduced densities, the correlation signal r](t) 
measured in || assumes the simple form 

v (t) =X { "\t)-X ( f\t). 

To the extent that the overlap r a (0) is and T a (t) remains negligible, it is clear from Eqs. 
(|i0|) that A_? ~ A+ , so that the above result reduces to 

r)(t) w \ { f\t) - X^it) = 1 - 2X i9 \t) w AL e) (t) - Af } (t) = 1 - 2X ( +\t) (11) 

which, together with the trace normalization, is sufficient to determine the eigenvalues in this 
small overlap limit. One sees moreover that, independently of the small overlap condition, 
the eigenvalues of the reduced densities are determined by the separately measured quantities 
P ee and P ge as A_ = P ee and A_ — P ge . 

Exactly the same analysis applies in this case when a master equation approximation 
is used to describe the damping effects associated with Eqs. ®, as done e.g. in ref. ||. 
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In fact, the master equation can be solved in closed form giving for the reduced density an 
expression identical in form to Eq. (|9|) with 



a (t) _> a (0)e~^ 2 and F b (t) - e -2|-(o)l 2 (i-e- <) 

where 7 = l/t c , the inverse of the damping time of the field intensity in C. 

Case b), which corresponds to the conditions of the experiment reported in ||, can be 
handled in precisely the same way as case a), since the state of the field prepared in C (i.e., 
immediately after the detection of the first atom) is still of the form of a superposition of 
just two coherent states, namely 

I t = 0} = N a(j) [e i(t> I ae^)+ | ae^)] <8> | b ). 

Expressions for the two non-vanishing eigenvalues and for the corresponding eigenvectors 
of the reduced densities can still be calculated in a straightforward way, in terms of the 
appropriate solutions of Eqs. (§), taking into account the non-orthogonality of the coherent 
state representation through the appropriate overlap matrix. Since they are considerably 
more cumbersome than in the previous case and will not be needed explicitly in the following 
development we do not give them in full here. 

Before turning to a discussion of the general case, we consider the simpler expressions 
appearing in the small overlap limit (i.e. | (a(t)e 1 ^ \ a(t)e~ l< ^) |<C 1 so that the overlap 
matrix essentially reduces to the unit matrix). Although the operators involved in U)jU e are 
no longer simultaneously diagonal with the reduced densities, one still gets for the correlation 
signal a result similar to Eq. (11), which can be written as 



V (t) « cos (^2 \Pk(t)\ 2 sin 2^ (A+ - a_; 



where A± ~ \± ~ A^ stand for the small overlap approximations to the eigenvalues of the 
reduced densities 

l±\T b (t,(/>)\ 
A ± — ^ 



with 

T b (t,cf>) = e -2Ej&WI 2 ^ 2 V e(t) and 6(t) = J2\Pk(t)\ 2 sm24>. 

h 

In the last equations a(t) and /3k(t) stand for the solutions to Eqs. (§) with initial conditions 
a(0), /3fe(0) = 0. This makes all the dependence on explicit. The possibility to express the 
small overlap approximation to r\ (t) in terms of reduced density eigenvalues is due to the fact 
that in this limit the two densities pf and p F g essentially commute and also have essentially 
the same eigenvalues, i.e. 
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pf «| A + )A + (A + I + I A_)A_(A_ | ; p F g «| A + )A_(A + | + | A_)A + (A_ | 

with approximate eigenstates 

| A±> oc ± e W®++) | a(t)e i<f> )+ | a(t)e"^). 

The correlation signal measured in || , Eq. @ , thus essentially determines the eigenvalues 
of the reduced densities pf and pf in the small overlap limit. Essential for this result is the 
pairwise equality of the eigenvalues of the two matrices in this limit. We then finally consider 
the general case and show that, independently of the small overlap assumption, separate 
measurements of P ee and P ge again determine the eigenvalues of the trace normalized matrices 
pf and pf . In fact, it is easy to show that these quantities provide in general independent 
linear combinations of the two non-vanishing eigenvalues in each case, with coefficients which 
depend on the dynamics of the atom-field interaction as expressed by matrix elements of the 
reduced operators We, Eq.(U), in the reduced density eigenvectors. In order to see this 
explicitly, it suffices to write the reduced densities in diagonal form as 

pf =| \f(t))\T(t)(\?(t) | + | £\t))\®{t)()P(t) | (12) 
which gives for the correlation functions 

P % % = (4\t) I U\U e | \?(t))\f(t) + (X-\t) I U\U e I \®(t))\®(t). 
These general expressions hold of course also in case a) above, where one has 

(Xf | U\U e | \f) = and (A® | U\U e \ = 1. 

As a result of the fact that in case a) the reduced densities are simultaneously diagonal with 
the reduced operators Ue , the determination of the eigenvalues becomes in fact independent 
of the structure of the (in this case, common) eigenvectors. 

We conclude from these results that the validity of the purely dispersive coupling as- 
sumption of the two-level atom to the field in cavity C, which leads to the validity of Eqs. 
(H|), and the adequacy of a coupling of the form of Eq. (|]) to implement the dissipative 
dynamics allows for a direct experimental monitoring of the intrinsic decoherence processes 
undergone by the field as a result of its coupling with the "heat bath". This can be done 
simply through the separate measurement of the conditional probabilities P ee and P ge , since 
they fully determine the non-vanishing eigenvalues (in this case only two, for all times after 
the preparation of the superposition state) of each one of the reduced densities pf and p g . 
Under less stringent conditions, allowing for measurable effects of higher order corrections 
(in fl/5) to Eqs. (||), one will generate field states, after the passage of the first atom, which 
cannot be reduced to superpositions of just two linearly independent states (see e.g. ref. 0]), 
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and will consequently have reduced densities having a correspondingly higher number of non- 
vanishing eigenvalues. The representation corresponding to Eq. (|T2| ) will then have more 
than just two terms so that the complete determination of intrinsic decoherence properties 
of the field will require fancier measurements. 
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